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ABSTRACT 

We investigate the evolution of the galaxy population since redshift 2 with a focus on the colour 
bimodality and mass density of the red sequence. We obtain precise and reliable photometric 
redshifts up to z = 2 by supplementing the optical survey COMBO-17 with observations in four 
near-infrared bands on 0.2 square degrees of the COMBO-17 A901-field. Our results are based 
on an ii-band-selected catalogue of 10692 galaxies complete to H = 21™7. We measure the rest- 
frame colour (C/280 — V) of each galaxy, which across the redshift range of our interest requires no 
extrapolation and is robust against moderate redshift errors by staying clear of the 4000A-break. 
We measure the colour-magnitude relation of the red sequence as a function of lookback time 
from the peak in a colour error-weighted histogram, and thus trace the galaxy bimodality out to 
z ~ 1.65. The (E/280 — V) of the red sequence is found to evolve almost linearly with lookback 
time. At high redshift, we find massive galaxies in both the red and the blue population. Red- 
sequence galaxies with log M* /Mq > 1 1 increase in mass density by a factor of ~ 4 from z ~ 2 to 
1 and remain nearly constant at z < 1. However, some galaxies as massive as logAf^/Af© = 11.5 
are already in place at z ~ 2. 



Subject headings: galaxies: evolution — galaxies: high-redshift 
stellar content — surveys 



galaxies: photometry — galaxies: 



1. INTRODUCTION 

The Lambda Cold Dark Matter (ACDM) cos- 
mological model predicts a bottom-up hierarchi- 
cal scenario of structure formation in which large 
structures have been formed recently (z < 1) 
through mergers of existing smaller structures. It 
is still unclear when and how in this context galax- 
ies have formed and how they have evolved into the 
baryonic structures made of stars we see today. 
The most evolved and massive systems known 
in the Universe are found especially among red- 
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sequence galaxies. They are straightforward to 
identify and seem to mark an end point in galaxy 
formation. Model predictions can thus be tested 
against measurements of the build-up of this popu- 
lation through cosmic time. To this end, the popu- 
lation needs to be tracked towards higher redshift. 
Their colour, mass and number density hold clues 
about their mass assembly over time. 

It is still unknown when in the history of 
the universe the red-sequence population first 
emerged. Using the optical COMBO-17 data, 
Bell et al.l l|2004 ) have shown that the galaxy bi- 
modality is present at all redshifts out to z = 1 and 
that the re d sequence was alrea dy in place at z = 
1, see also IWeiner et al. I (|2005h . Several spectro- 
scopic stud ies have observed massive red galaxies 
at z > 1.5 dCimatti et al. 120041 Idazebrook et al 
2004; iDaddi et al.l 120051 iMcGrath et al.1 12007: 
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Conselice et al.ll2007 ; Kriek et al . 2008: Cass ata et al 
20081 ). but did not allow the study of large sam- 
ples. On the other hand, photometric redshift sur- 
veys combining optical and NIR data can provide 
photometric redshifts for lar ger samples o f high 
redshift galaxies. Recently, Taylor et al. (2009) 
demonstrated the persistence of th e red s equence 
up to z ~ 1.3 and ICirasuolo et al. ( 2007 ) as well 



as iFranzetti et al.l (|2007l ) up to z ~ 1.5. Us- 
ing colour-colour diagrams ([Williams et al. 2009) 
claimed the presence of quiescent galaxies anal- 
ogous to a red sequence up to z ~ 2. A per- 
sistent issue for photometric redshift surveys is 
an increasing uncertainty in the measurement 
of rest-frame colours at high redshift (z > 1.5). 
It increases the scatter in colour-magnitude and 
colour-stellar mass diagrams, thus blurring the 
signature of a possible red sequence at high red- 
shift. Especially, the expected small red-sequence 
population in field samples may be undetectable. 

In this paper, we maximise the separation of red 
quiescent and blue star-forming galaxies by using 
the rest-frame co lour (E/28O — V) as a, diagnostic. 
Wolf et al.l (|2009i) have shown how in galaxy clus- 
ters this colour also tends to render star-forming 
red galaxies bluer to leave a more purified red sam- 
ple. We use colour-error-weighted histograms to 
improve the recognition of a red sequence in the 
face of increasing colour uncertainties at high red- 
shift (see Sect. 3 for all methods). In Sect. 4 we 
investigate the colour distribution up to z = 2 and 
derive a colour-magnitude relation and a bimodal- 
ity separation evolving with redshift. Thus, we 
attempt to estimate when the red-sequence popu- 
lation emerges. Ultimately, we quantify the evo- 
lution of the following galaxy properties: colour, 
luminosity, mass and number density for the red- 
sequence galaxy population since z = 2 (Sect. 5). 
We summarize and conclude this study in Sect. 6 
Throughout this paper, we use the cosmological 
parameters fl m = 0.3, Oa = 0.7 and H — 70.7 
km s" 1 Mpc" 1 . All magnitudes are in the Vega 
system. 

2. OBSERVATIONS 

2.1. The COMBO-17+4 Survey 

The COMBO-17+4 survey is the NIR exten- 
sion of the optical COMBO-17 survey and is de- 
signed to probe galaxy evolution since z = 2. Near 



infrared data were necessary to obtain optical rest- 
frame properties for galaxies in the redshift range 
1 < z < 2, where the optical rest- frame of galax- 
ies is shifted into the NIR. The survey consists 
of observations in four NIR bands (A/AA): the 
three medium-bands Y (1034/80), J x (1190/130), 
J 2 (1320/130), and the broad-band (1650/300). 
In the long run, the COMBO-17+4 survey targets 
three independent fields (A901, A226, and Sll) 
for a total coverage of 0.7D° (see Table [T] for co- 
ordinates and integration times). The results de- 
rived in this paper are based on the observations of 
the A901-field only. This field contains the super- 
cluster of galaxies A901/2 located at z = 0.165 
where further multi-wavelength coverage from X- 
ray to radio was o btained by the STAGES survey 
dGrav et al.ll2009h . 



2.2. Data 

The NIR data were obtained in several ob- 
serving runs from December 2005 to April 2009 
with the NIR wide field camera Omega2000 at the 
prime focus of the 3.5-m telescope at Calar Alto 
Observatory in Spain. The camera has a pixel size 
of 0.45" and a wide field of view of 15.4' x 15.4', 
so that a half-degree COMBO-17 field can be cov- 
ered with a 2 x 2-mosaic. On the A901-field, only 
three of four pointings could be finished in the time 
awarded to the project, and hence the NIR data 
are missing in the south-west quadrant. 

Also, an area of 68" x 76" centred at (a, S) J2000 = 
09 ?i 56 m 32 s .4, -10°01'15" was cut out to avoid 
spurious objects created by the halo of a very 
bright (K = 5™75) Mira variable star, also known 
as the IRAS point source 09540-0946. The to- 
tal NIR coverage is thus 690D' ~ 0.19D . The 
optical data include a combination of 17 broad 
and medium bands centred between 365 nm and 
915 nm and were obtained with the Wide Field 
Imager (WFI) at La Silla Observatory between 
February 1999 and Jan uary 2001 by the C OMBO- 
17 survey (|Wolf et alJl2003t foray et alJ l2009). 



2.3. Data Analysis 

The NIR data reduction and photometry were 
performed using the software ESO-MIDAS in com- 
bination with the MPIAPHOT pack age devel- 
oped by Ro ser fc Meisenheimer ( 1991 ) and the 
OMEGA2k data reduction pipeline developed by 
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Fafibenderl ([20031 ) . The image data reduction con- 
sisted of flatfielding, dark current and sky back- 
ground subtraction as well as correction for bad 
pixels and pixel hits by cosmic ray events. Ad- 
ditive stray light appearing in a ring shape in all 
scientific and calibration images taken with the 
Omega2000 camera have been subtracted at the 
flatfield level for images taken in the Y, J\ , and J 2 
filter where the stray light, if not corrected, would 
have contributed an additive 5%, 5%, and 10% to 
the flux, respectively. No stray light correction 
has been done for the iJ-band images where the 
additive contribution is negligible at < 0.5%. 

From the coadded £f-band images of the three 
pointings we created an ii-band mosaic image 
with a total exposure time of 11600 sec pixel -1 
and 1.03" seeing on average. The summation pro- 
cess assigned a weight to each input image accord- 
ing to its transmission, background noise and PSF, 
so that the g-ban d mosaic has an optimal PSF. 
Using SExtractor (jBertin fc Arnoutsl Il996l ) with 
default parameters we obtained a deep i7-band 
source catalogue with 31747 obje cts. The ast rom- 
etry was performed with IRAF (jTodyl 119931 ) us- 
ing hundreds of bright (H < 16™0) 2MASS point 
sources in common with our catalogue reaching an 
accuracy of 0.1" in RA and DEC. 

The optical photometry has been re-derived for 
the source positions in the -ff-band catalogue in- 
stead of the previous i?-band catalogue. Hence, 
optical and NIR photometry in all 21 bands of 
the COMBO-17+4 survey are measured in aper- 
tures matched in location and in size. In prac- 
tice, we use a Gaussian weighting function in our 
aperture to give more weight to the bright central 
parts of an object and less weight to the fainter 
outer parts. Using a Gaussian sampling function 
and assuming a Gaussian seeing PSF allows us to 
sample identical areas of an object independent of 
seeing: we adjust the width of the Gaussian aper- 
ture to counteract seeing changes such that the 
integral over the aperture is conserved under see- 
ing changes. Mathematically, our brightness mea- 
surements are identical to placing apertures with a 
Gaussian weighting function of 1"7 FWHM onto a 
seeing-free image in all bands (COMBO-17 chose 
1"5 given the slightly better seeing of its optical 
data). The NIR filters Y, J x , J 2 and H reach 10- 
cr aperture magnitude limits of 22™1, 21™5, 21^4, 
and 21'. n 0, respectively. 



The COMBO-17 spectrophotometry stars have 
been used to ensure the c alibration of the op tical 
bands between each other (jWolf et al.ll2001bl ). We 
extended the calibration into the NIR using the 
Pickles! ([19981 ) spectral library by visually match- 
ing the colours of point sources in our data to main 
sequence stars in the library. We estimate that the 
relative calibration between optical and NIR has a 
limited accuracy on the order of 7%. The calibra- 
tion of our photometry in the H-, J 2 -, and Ji-band 
has been verified by comparing the H magnitude 
and the average value of the J2 and J\ magnitude 
respectively to the H and J magnitudes of 340 
stars with H < 16 in common with the 2MASS 
catalogue. For the point sources, we found a mean 
offset in magnitude of < 1% in the H-b'&nd and 
< 0.4% in the averaged J\- and J2-band. No ver- 
ification could be done for the F-band magnitude 
since no point source catalogue currently exists in 
this waveband. 

3. METHODS 

3.1. Photometric Classification and Red- 
shifts 

Photometric redshifts were determine d us- 
i ng th e multi-colour classification code by IWolfl 
(| 19981 ) as it was done in COMBO-17. Objects 
are divided into the four classes star, white dwarf, 
galaxy and quasar by comparing measured colours 
with colour libraries calculated from spectral tem- 
plates. The template s for stars and qu asars are 



identical to those in Wolf et al.l ( 2004 ). while a 



modified library has been built for galaxies, in or- 
der to allow reliable estimates of the stellar mass 
(see below). At the bright end, the redshift accu- 
racy is limited by systematic errors in the relative 
calibration of the different wavebands or in a mis- 
match between templates and observed spectra. 
At the faint end, photon noise dominates. 

This dataset is a superset of the COMBO-17 
data with the four NIR bands added. As a result, 
the photo-z's have changed very little at z < 1, 
where the NIR bands add little constraints, but it 
is reasonable to assume that they have improved 
over the optical-only results at z > 1. In Fig. [T] 
we show four examples of red galaxy SEDs and 
their best-fit templates. At z < 1 the fit is clearly 
constrained by the original COMBO-17 data at 
A < 1/im. At z > 1 the four NIR bands act in 



3 



concert with the optical data, while for red galax- 
ies towards z = 2 they are the sole providers of 
significant flux detections. The two z > 1-galaxies 
shown are examples of EROs with R — H ~ 5 
and 5.7, respectively. To the eye, their redshift is 
clearly constrained by locating the break between 
neighbouring pairs of filters, while the fit takes all 
filters into account to constrain the redshift fur- 
ther. 

COMBO-17 photo-z's have been shown to be 
accurate to aj(l + z) < 0. 01 at R < 21, < 0.02 
at R < 23 (jWolf et alJl2004h . albeit on a different 
field (the CDFS). On the A901 and Sll field, we 
only have spectra for galaxies at z < 0.3. From 
these, the photo-z dispersion of t he cluster A901/2 
has been measure d as 0.005 rms ( Wolf et al]|2005t 
Gray et al.ll2009l) . Presently, we lack the ability 
to confirm our photo-z's at z > 1 here. Hence, we 
need to rely on the plausibility of the SED fits to 
the photometry as shown in Fig. [TJ 

However, the photometric errors and the grid of 
templates allow the estimation of probability dis- 
tributions and confidence intervals in redshift for 
each galaxy. These estimated redshift errors are 
shown in Fig. [3J where they illustrate the change 
in behaviour with redshift and magnitude. Ob- 
viously, photometric errors increase towards faint 
magnitudes and propagate into redshift errors, but 
at z < 1.2 the redshift error is still very much 
driven by the deep optical photometry. Hence, 
the bulk of objects has er 2 /(l + z) < 0.05 even at 
our adopted ff-band limit. At z > 1.2, however, 
where the main redshift constraints are in the NIR 
bands, a clear upswing of redshift errors with H- 
band magnitude can be seen. 

3.2. Galaxy Samples 

After eliminating ~ 2300 objects classified as 
non-galaxies, we defined our sample by first apply- 
ing an ii-band selection above our 5-er detection 
limit of H = 21™ 7. However, we are concerned 
with the evolution of galaxy samples that need to 
be complete in rest-frame l/-band luminosity or 
stellar masses. We avoid a colour bias by apply- 
ing further magnitude cuts and eliminate from our 
colour-magnitude diagrams those faint tails of the 
galaxy distribution that are known to be incom- 
plete. Across our redshift range several observed 
bands with different completeness limits map onto 
the rest-frame y-band. Hence, we opted to apply 



the following two further cuts: 

(1) At z < 0.43 the entire observed i?-band 
is redwards of the 4000A-break and roughly co- 
incides with the rest-frame y-band. Here, we re- 
quire galaxies to have R < 23 m 5, which is the 
completeness limit of the optical-only COMBO- 
17 redshifts for z < 0.43-galaxies. The presence of 
NIR data may have deepened our completeness in 
this regime, but we opt to err on the conservative 
side and eliminate this faint end from our analysis. 
Note, that optically-faint red galaxies at higher 
redshift (such as EROs) are by definition NIR- 
bright and have thus very well-constrained SED 
fits and (higher) redshift estimates (see Fig. Q}. 

(2) In the redshift range 0.43 < z < 1.4, the 
Y"-band is redwards of the break and we require 
Y < 22 m 8 (our 5-cr-limit). Again, objects that are 
particularly red in Y — H are expected to reside 
at z > 1.4. 

Finally, at z > 1.4 the SED around the 4000A- 
break is entirely sampled by our NIR filters and 
the H < 21™ 7 selection is sufficient to ensure 
completeness and an accurate redshift estimate. 
These selections are complete in stellar mass 
log(M*/M ) to 8.5, 9.5 and 10 in the redshift 
ranges z < 0.43, 0.43 < z < 1.4, and z > 1.4, 
respectively. 

After removing galaxies with bad flags, our 
sample contains 10692 galaxies. Fig. [3] shows this 
sample and reveals already a number of large-scale 
structures. The dominant overdensity at z ~ 0.16 
is the original main target of this COMBO-17 
field, the supercluster region A901/2 with well- 
over 1000 galaxies. Further clusters and large- 
scale structures have been identified and reported 
i n this field, partly f rom weak gravitati onal lensing 
(jTavlor et al.ll2004 Isimon et alj|20ld) and partl y 
from a galaxy cluster search ( Falter et al. 2010h . 
These include localised clusters embedded in large- 
scale structures at z ~ 0.26, 0.37, 0.5, 0.7, 0.8. At 
higher redshift and fainter magnitudes, i.e. z > 
0.8 and H > 20, redshift focussing effects are pos- 
sible, and structures which appear only at faint 
magnitudes but without a sharp tail to the bright 
edge are unlikely to be physically localised over- 
densities. In Fig. [3] the two possible unreal over- 
densities are the faint-end blobs at z ~ 1 and 
z ~ 1.2. Their appearance does not require red- 
shift errors in excess of what is discussed above, 
but only a mild focussing within the allowed er- 
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rors. At z > 1.2 no structure can be seen, as any 
physical contrast has been smoothed by our red- 
shift errors. 

The histogram of photometric redshifts of the 
complete galaxy sample is shown in Fig. [4] Here, 
we show first the optical-only COMBO-17 red- 
shifts of four different fields, which demonstrate 
the signature of abundant large-scale structure. In 
case of the CDFS, these have been reported con- 
sistently by a variety of groups. Overall, the struc- 
tures include six Abell clusters as well as further 
clusters and rich groups, usually characterised by a 
pronounced red sequence. It is thus clear that we 
investigate a variety of environments across red- 
shift, and only the combination of several fields 
will eventually suppress error propagation from 
ficld-to-field variation. 

The histogram in the right panel shows new red- 
shifts after including the NIR bands of COMBO- 
17+4. It is selected by combined apparent- 
magnitude cuts in the R-, Y- and _ff-bands, in- 
stead of the sole i?-band in the optical-only photo- 
z's, and thus numbers at higher redshifts are in- 
creased. The inclusion of NIR data allows for 
good photo-z's at z > 1, but presently such data 
are only available for three quarters of the A901- 
field covering 0.2D° area. This is our first look at 
z > 1.2-galaxies, while completion of our project 
will entail three fields with a different mix of en- 
vironment at each redshift. 

The bulk of our sample is in the range 0.7 < 
z < 1.1, consistent with the z = 0.9-peak in the 
n(z) found in COSMOS (jilbert et al.ll2010l) . 



on a two burst-mo del following the approach of 



3.3. Spectral Library of Galaxies 

The photometric redshi fts of galaxies in the 
optical-only COMBO-17 (|Wolf et al.l 120041 ) are 
based on synthetic templates prod uced by the 



Borch et al.l (|2006l ). The motivation for a two- 
burst model is to reproduce better the stellar 
population of blue galaxies that contain both an 
old stellar population and ongoing star-formation 
activity and to provide more realistic mass-to- 
light ratios. How ever, the two-burst model by 
Borch et al. ( 20061 ) failed to deliver accurate red- 



popu lation synthesis model PEGASE (jFioc fc Rocca 
19971 ). A single, exponentially declining burst of 
star formation was assumed for the star forma- 
tion history (r = lGyr). The templates form a 
2-dimensional grid in which one parameter is the 
age since the start of the burst (60 steps between 
0.05 to 15 Gyr) and the second one is extinction, 
which was modelled as a foreground screen of dust 
following the S mall M agellanic Cloud (SMC) law 
defined bv IPeil ( 19921 ) in six steps from Ay = to 
Ay — 1.5. Here, we employ a new library based 



shifts. 

Our new two-burst library solves this problem 
by including extinction again and adjusting the 
age and strength of the second burst. For red 
galaxies (ages > 3 Gyr) we leave the single-burst 
templates unchanged (t\ = lGyr), and for bluer 
galaxies we add a recent burst (ji = 0.2 Gyr) to 
the old population starting 2.75 Gyr after the first 
one. Increasingly blue galaxy templates are gener- 
ated by both increasing the relative strength of the 
second burst and by moving the final age (at which 
the galaxy is observed) from 3.0 Gyr to 2.80 Gyr 
(i.e. towards the start of the second burst). This 
ensures that the templates fall into a region of the 
rest-frame colour-colour diagrams that is occupied 
by observed galaxies, and thus leads to very accu- 
rate photometric redshifts. We assume a Kroupa 
IMF, an initial metallicity of 0.01 (~ 2/3 Z Q ) and 
neither infall nor outflow. More details will be 
given in Meisenheimer et al. (in preparation). 

3.4. Stellar Mass Estimation 

The stellar mass of a galaxy is derived from 
the best-fitting template SED, whose characteris- 
tics are constrained mostly by the rest-frame spec- 
trum between 280 nm and the V^-band, which is 
observed at all redshifts. Formally, we use data 
from all bands, but our masses can be seen to 
be derived from a galaxy's observed y-band lu- 
minosity and the 1^-band mass-to-light ratio of its 

-VolmMeT ^g tem P late - 



Our use of a fixed rest-frame band to deter- 
mine a gala x y's lu minosity differs from that of 
Borch et al. (|2006l) who only used the reddest 



observed-frame band with good photometry, and 
thus may suffer unintended bi as trends with red - 
shift. At low redshifts z < 0.5. lBorchet al.l (|2006l) 
sample the rest-frame V-band as well so that our 
mass estimate should only differ from theirs due 
to the different templates. Only due to the NIR 
extension in COMBO-17+4 can we continue to 
use the rest-frame F-band out to z = 2. 
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Indeed, a comparison of galaxies at the red- 
shift of the super-cluster Abell 901/2 (0.150 < 
z < 0.175) shows that systematic differences in 
the stellar masses do not exceed 0.1 dex. Beyond 
z = 0.5 the stellar masses derived here are supe- 
rior due to the NIR photometry. We estimate a 
typical mass accuracy on the order of 30% across 
all redshifts. 

We wish to assess possible biases in the stellar 
mass estimates arising from dust and consider the 
Bell fc de Jong I ([20011 ) equation, 



log 10 M*/M oc -0AM v + 1.305(5 - V) (1) 

because it allows an analytic derivation of a red- 
dening vector in colour-mass diagrams, and our 
masses are consistent with Eq. 1 within ±0.1 dex. 
Any dust reddening can now be seen as a mean 
reddening plus some structure. Approximating 
the mean reddening by a uniform screen of dust 
right in front of the stellar population changes its 
colour by A(B — V) = Eb-v and its F-band 
luminosity by AV — RvEb-v- The reddening- 
induced overestimate of the M/L-ratio and the 
absorption-induced underestimate of L cancel ex- 
actly, if 



0ARvEb-v = 1.3055s- 



Rv = 3.2625 (2) 



which is almost tr ue for the Milky Way, LMC 
or SMC dust law from|Pe3 (|l992l ). Dust acting like 
a foreground screen on the stellar population will 
thus not bias the mass estimates at all. However, if 
pockets of highly absorbed stars exist as well, they 
will be entirely withdrawn from the optical view 
and not contribute to either luminosity or colour. 
They will thus be plainly not taken into account 
and the final mass estimate will be underestimated 
by just the stellar mass present in highly obscured 
pockets. 

3.5. Rest-Frame Luminosities and Colours 

Rest-frame luminosities are derived from the 
observed photometry covering the wavelength 
range from the t/-band to the iJ-band. We ob- 
tained rest-frame luminosities for the 1^-band in 
the Johnson photometric system as well as for the 
^280-band, a synthetic UV continuum band with 



a top-hat transmission curve that is centred at 
A = 280 nm and 40 nm wide. These two filters are 
covered by our observed bands across almost the 
entire range of interest. This is essential as our 
study does not rely on SED extrapolation except 
for galaxies in the small range 0.2 < z < 0.3. 
The rest- frame colour [/ 2 so — V is measured ro- 
bustly against redshift errors since both filters 
are located in smooth continuum regions of the 
galaxy spectra. The Johnson [/-band, in contrast, 
partly overlaps with the 4000A-break and is hence 
affected strongly by small uncertainties in the red- 
shift determination. 

The rest -frame lumi n ositie s are derived as de- 
scribed by IWolf et all ([20041 ). The best-fitting 
SED is placed into the aperture photometry and 
integrated over the redshifted rest-frame bands. 
We have taken into account the interstellar fore- 
ground extinction as well as a correction from 
aperture to total photometry, which could be bi- 
ased by colour gradients. It is determined from 
the total magnitude MAG-BEST derived by SEx- 
tractor on the ii-band image and thus certainly 
correct for any rest-frame band overlapping with 
the observed -ff-band. The magnitude errors are 
determined from a propagation of photometric er- 
rors. They include a minimum error of 0™1 to take 
into account redshift errors and overall calibration 



uncertainties, for details see IWolf et al.l (J2004J). 



3.6. Error- Weighted Colour Histograms 

Generally, colour and magnitude information 
can be used to separate passive red galaxies from 
star forming blue ones. However, at high red- 
shift the well-known colour bimodality may be 
blurred by scatter from relatively large errors on 
the colour measurements. We try to manage 
this challenge and investigate the galaxy colour 
bimodality through cosmic time by using error- 
weighted colour histograms. This method repre- 
sents each galaxy by its Gaussian probability dis- 
tribution in colour c, i.e. p(c) ~ e - ^ 00 -* / 2cr ^, 
where a c is the colour error and p is normalized 
so that J pdc=l. Thus, a galaxy with a small 
error has a more peaked distribution and con- 
tributes more structure to the summed distribu- 
tion than a galaxy with a large error. As a re- 
sult, the structure in our histograms is driven en- 
tirely by galaxies with small colour errors, while 
objects with large errors lift the overall counts 
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without producing peaks. We produced error- 
weighted colour histograms by summing all the 
Gaussian distributions within many thin redshift 
slices (Az ~ 0.1) stepping through our full redshift 
range of < z < 2. Towards highest redshifts, the 
increasing colour errors will dilute the contrast of 
the red-sequence peak. Redshift errors can also 
lead to some spill-over from a physical redshift bin 
into neighbouring photo-z bins and produce scat- 
ter in our measured colour evolution from bin to 
bin in redshift. 

To better disentangle the red sequence from the 
blue cloud we have tilted the (U280 — V) colour 
in the colour-magnitude plane, see Eq. [3l The 
measured (L^so — V) of each individual galaxy is 
projected along the slope of the red sequence as de- 
termined in the CMD of the super-cluster A901/2 
(see Fig. [5J top left) to the pivotal magnitude 
M v = -20: 



((7 28 o - V) Uv =-20 = ((7 28 o - V) + 0.3(M V + 20) (3) 



The derived slope of 0.3 is co nsistent with tha t 
of 0.08 in the (U - V) colour in lBell et all (|2004f ). 



4.1. 



EVOLUTION OF THE COLOUR BI- 
MODALITY 

Rest-Frame Colours and Colour Er- 
rors as a Function of Redshift 



Scatter in measured galaxy colours affects the 
appearance of the colour bimodality in a colour- 
magnitude diagram and may render it invisible, 
especially in a sample of low density field galax- 
ies where no clusters rich in red galaxies produce 
a well-defined sequence. As the scatter in colour 
increases with redshift, this effect can prevent us 
from detecting a physically present high-redshift 
bimodality. 

For a first assessment of our ability to trace bi- 
modality, we plot the rest-frame (L^so — V) colour 
and its error 5(U2so ~ V) as a function of redshift 
in Fig. [5] In panel (a), we see the number of red 
galaxies decreasing with redshift. Also, a lack of 
star-forming blue galaxies with (?/280 — V) < 1 at 
z > 1.5 is caused by a red-object bias from our 
ii-band selected catalogue. Panel (b) shows that 
the colour accuracy is very good for objects lo- 
cated at low redshift z < 0.9. The bulk of the 



galaxies have a colour error close to the assumed 
minimum of 0™1. The large scatter around z ~ 1 
is caused by a locally increased magnitude error in 
the rest-frame V-band that results from My be- 
ing calculated from the relatively shallow narrow- 
band J\. In contrast, the C/280 filter overlaps with 
the ~ 3 mag deeper i?-band at z ~ 1. At z > 1.7, 
the colour uncertainties grow larger again since 
the rest-frame where the L^so-band falls into the 
/-band, which is our shallowest broad-band. 

4.2. Colour Bimodality and the Emer- 
gence of the Red Sequence 

In Figure [5] we plot error- weighted colour his- 
tograms of galaxies in different redshift slices. 
Each panel shows two distinct peaks due to the 
bimodality, whereby the right peak represents red- 
sequence galaxies sitting on top of a tail of blue 
cloud galaxies reaching smoothly towards the red 
due to lower star-formation rates or dust red- 
dening. Clearly, a red galaxy sample defined by 
a colour cut will comprise both quiescent red- 
sequence galaxies as well as dusty star-forming 
galaxies. At high redshift, we are unable to sep- 
arate old red from dusty red galaxies, and so a 
red galaxy sample overestimates the space den- 
sity of quiescent galaxies. However, the colour- 
magnitude relation should be well constrained, as 
only galaxies from the proper quiescent red se- 
quence are focussed in colour and contribute to 
the peak in the colour histogram, while dusty 
red galaxies form a smooth underlying continuum 
spreading across the red-sequence and extending 
beyond. 

The super-cluster A901/2 seen in the top left 
panel has a particularly clear red sequence. Our 
results show a clear gala xy bimoda l ity at z < 1 
as already established by iBell et al.1 (l2004h . How- 
ever, the depth of our survey allows us to ex- 
tend the detection of the galaxy bimodality up to 
z ~ 1.65, the mean redshift of the highest inter- 
val where we can still detect two distinct distri- 
bution peaks. Beyond z ~ 1.65 it is not possible 
to confirm a galaxy colour bimodality since the 
number of objects available for our analysis drops 
considerably. Figure[6]bottom right shows the dis- 
tribution of 423 galaxies in the redshift interval 
1.7 < z < 2. It is clearly not bimodal, though 
this does not mean that the red sequence does not 
exist there. Our redshift errors are unlikely to be 
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the main source of this disappearance, since the 
redshift error cr z /(l + z) of luminous red galaxies 
with My = [—23.5, —22] degrades only from a me- 
dian value of ~ 0.01 at z = 0.8 over ~ 0.025 at 
z = 1.2 to ~ 0.04 at z = 1.6. Instead, it is the 
combination of a relatively small number of qui- 
escent luminous red objects due to the small area 
surveyed and the washing out of any red sequence 
signal due to increasing colour errors (see Fig. [SJd) 
that prevents us from detecting the red sequence 
beyond this redshift. 

4.3. Evolution of the Colour-Magnitude 
Relation 

We use the error-weighted colour histograms 
in thin redshift slices of Az ~ 0.1 to track the 
colour evolution of the peak in the red sequence 
and plot (U 2 so — V) at the pivot point M v — —20 
in Figure H The colour of the red-sequence peak is 
obtained with the MIDAS routine CENTER/GAUSS 
gcursor, which fits a Gaussian function to an 
emission line on top of a continuum. We mark 
the 'continuum' level on both sides of the peak 
interactively, and obtain the peak location from 
the fit. We obtain errors from the uncertainty of 
the peak position obtained by varying the inter- 
val in which the Gaussian is being fit, but we do 
not include systematic calibration uncertainties in 
colour (see Table [5] for results). We find that the 
(^280 — V)m v =-20 colour of the peak evolves al- 
most linearly with lookback time t, and obtain a 
linear best fit (dashed line) of 



In the low redshift regime < z < 1 this cut 
i s con sistent with the one derived by Bell et al 



280 



V) Mv =-20 = 2.57 - 0.195 x r (4) 



Our data points of the peak colour scatter 
somewhat around this fit, which is an indication 
of our uncertainties in measuring the colour, as 
any cosmological evolution of the average galaxy 
population is expected to be monotonic. However, 
these variations may also result from us looking at 
different environments at different redshift, which 
may be at different evolutionary stages, and not 
solely from our methodical uncertainties. 

We obtain a colour-magnitude relation (CMR) 
as a function of redshift using the approximation 
for the lookback time r ~ \hGyr ■ z ■ (1 + 2) -1 , 
and separate blue and red galaxies with a parallel 
relation 0.47 mag bluer than the CMR, which is 

{U 2 m-V) Um = 2.10-0.3(My+20)-2.92z/(l+z) (5) 



(J2004J) considering an approximate colour trans- 
formation derived from our templates, which is 
(U - V) = 0.28 + 0.43([/ 280 - V). Our CMR is 



slightl y steeper than the one derived bv lBell et al 
( 20041 ). but that translates only into a minor 



colour difference of < 0.1 mag for galaxies in the 
relevant magnitude range — 24 < My < — 18. Al- 
though we see no clear red sequence at z > 1.65, 
we extrapolate the CMR all the way up to z = 2 
to isolate red galaxies across our entire sample. 

The solid line in Figure [7] shows the evolution 
predicted by PEGASE for a single stellar popu- 
lation formed at a lookback time of 12 Gyr (zf 
= 3.7) with solar metallicity. The bulk of our 
data are consistent with the pure aging model. 
Our lowest-redshift value is based entirely on the 
A901/2 super-cluster at z = 0.165 (Lookback 
time=2.0 Gyr), which is the highest-density en- 
vironment in our sample. 

5. EVOLUTION OF THE RED-SEQUENCE 
GALAXY POPULATION 

We now focus on the colour and the mass evo- 
lution of the whole red galaxy population, and 
analyse colour-magnitude diagrams (CMD) and 
colour-stellar- mass diagrams (CM*D) in redshift 
slices. We divide our galaxy sample into red and 
blue populations with the cut in Eq. [5] Across 
the entire sample at < z < 2, roughly a third 
of the galaxies are red (3163 out of 10692), but in 
the high-redshift part at z > 1 less than a quarter 
are red (843 out of 3479 galaxies). 

5.1. Colour-Magnitude Diagrams 

The evolution of the red-sequence population is 
presented in the CMDs of Fig. [8] In each panel, 
the solid line indicates the bimodality separation 
of Eq. [5] assuming the mean redshift of the slice. 
The colour of the data points shows their individ- 
ual nature, and due to the width of the redshift 
intervals some galaxies scatter across the sepa- 
rating line. The top left panel in Fig. |8] shows 
the CMD of the A901/2 super-cluster centred at 
z = 0.165. Such a dense environment shows a 
clear red sequence, and this was used to derive 
the red sequence slope in Eq. [3] The red sequence 
is not as sharp and easy to visually distinguish 
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on CMDs beyond z = 1, especially in the redshift 
slice 1.07 < z < 1.19 due to the local increase 
in the scatter caused by large colour uncertain- 
ties. This highlights the necessity of using error- 
weighted histograms to derive the CMR. 

Altogether, we find that at a given magnitude 
both galaxy populations were bluer in the past, 
and in particular the bright end of the red se- 
quence became ~ 0.4 mag redder from z — 2 to 
z = 0.2. We also find an increasing population 
of bright (My < -22) blue ((U 280 - V) < 1) 
ga laxies at z > 1, wh ich has also been reported 
bv lTavlor et aD (l2009h in (U - R)-vs.-M R CMDs 
of a i^-selected sample at 0.2 < z < 1.8 in the 
Extended Chandra Deep Field South. 

5.2. Colour- Stellar Mass Diagrams 

Fig. [9] shows the evolution of the population in 
CM*D for the same redshift slices as the CMDs in 
Fig. [8l We select the red sequence (black points) 
again with Eq. [5] We reach smaller masses than 
Borch et"aH ( 2006 ) because our galaxy sample is 
primarily i7-band selected instead of i?-band se- 
lected. Conversely, our blue galaxy populatio n 
reaches less deep than that of iBorch et al.1 (|2006h . 
We see a general trend whereby at fixed mass both 
galaxy populations were bluer in the past, just as 
they were at fixed magnitude. We find that mas- 
sive galaxies between 0.2 < z < 1.0 are domi- 
nated by the red population . This was also ob- 
served by IBorch et al. ( 20061) . who derived mass 
functions for the red and the blue galaxy popula- 
tion at z < 1 using the optical CO MBO- 17 data, 
and b y many other authors (see e.g. iBundv et al. 
20051) . However, at z > 1 there is a growing 



population of massive (log M*/M Q > 11) blue 
(([^280 — V) < 1) galaxies as we g o back in time . 
a ph e nomenon (also ob served by IWilliams et al 



120091 : iTavlor et ai1l2009l ) that has no analogue in 
the local Universe. 

The bottom right panel shows our highest red- 
shift slice (1.78 < z < 2) where our sample con- 
tains 45 red galaxies. Among them are eight 
very massive (logM*/M Q > 11.5) objects, which 
means that at z ~ 2 the red sequence already 
contains very massive galaxies. The most mas- 
sive object found in this high redshift slice has 
a stellar mass of logM*/M = 12.0. Our red se- 
quence will contain dusty star-forming galaxies be- 
sides old galaxies, but since these do not form a 



sequence and are not focussed in colour, they will 
not affect the measurement of the red-sequence 
colour with our tailored method. 

5.3. Number Density Evolution of Mas- 
sive Red Galaxies 

Figure [10] shows the evolution of the number 
density of massive red galaxies as a function of 
redshift (see also Table [3]) . These red galaxies in- 
clude both quiescent old galaxies as well as dust- 
reddened galaxies. We avoid the super-cluster 
A901/2 by constraining our sample to 0.2 < z < 2 
and choose masses of log M* /M Q > 1 1 , where we 
are complete at all redshifts, thus retaining only 
478 objects. 

We find that the number density of massive 
red galaxies rises considerably (by a factor ~ 4) 
from z ~ 2 to z ~ 1 and is more or less constant 
at z < 1. Due to the small survey area our re- 
sults are affected by cosmic variance on the orde r 
of 30% as estimated using Moster et al. ( 2010l ). 
Nevertheless, our results are comparable within 
the error bars to results by the GOODS-MU SIC 
survey (0.04 deg 2 area) (jFontana et al.l 120061) and 



the M USYC survey (0.6 deg 2 area) (|Tavlor et al 
2009). In the highest redshift bin at z ~ 2 the 



number density der ived is consis t ent w ith the spec- 
troscopic survey of iKriek et al.l (|2008l ). Since any 
sample of red galaxies could in principle be con- 
taminated by some dusty red galaxies the num- 
ber density derived here represents an upper limit. 
However, we do not expect the effect to be large 
at the high-mass end. 

5.4. Stellar Mass Density Evolution of 
Massive Galaxies 

The different evolution of stellar mass density in 
red and blue galaxies is shown in Fig. [TT] (see also 
Table [3]) . Again we restrict ourselves to masses of 
log M*/Mq > 11, where we are complete across 
0.2 < z < 2. Like the number density, the stellar 
mass density of the red population increases by a 
factor of ~ 4 from z ~ 2 to z ~ 1 and remains 
roughly constant at lower redshifts. In contrast, 
the stellar mass density of the blue population 
increases only by a factor of ~ 2 from z ~ 2 to 
z ~ 1.2 and decreases by the same factor towards 
low redshift again. Altogether, the overall stel- 
lar mass density increases by a factor of ^3 from 
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z ~ 2 to z ~ 1 due to the combined contribution of 
red and blue galaxies, while it remains constant at 
lower redshift where the red sequence dominates. 
Thus, the main formation epoch of the massive 
red galaxy population is ranges over 2 > z > 1. 

Comparing our stellar mass densities with the 
literature is not foolproof since different surveys 
have different selection criteria using e.g. colour, 
morphological types and star formation activity. 
They employ different methods to derive the stel- 
lar mass and they can be affected by cosmic vari- 
ance. 

Our results for the evolution of the mass den- 
sity of the entire po pulation between 1 > z > 
are consistent with IConselice et al.l (|2007l ) who 



also found little evolution for galaxies with 11 < 
log M* /M@ < 11.5. At higher redshift, our results 
differ: IConselice et al.1 (|2007l ) found an increases 
by a factor 10.7 from z ~ 2 to z ~ 1, while our 
results show a more modest rise by a factor ~ 3. 

Based on a 3.6M- selected sample of galaxies 
Arnouts et al.l ( 2007 ) found an increase of the 
mass density in the quiescent population by a fac- 
tor of 2 from z ~ 1.2 to z ~ 0, while the star- 
forming population shows no evolution. Addition- 
ally, at higher redshift between 2 > z > 1.2 they 
found that the quiescent population increases by 
a factor of 10 while the star formation population 
increases by a factor of 2.5. However, these results 
are based on a magnitude-selected sample and not 
a mass-selected one. 

found in a sample of 
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Mk < — 23 galaxies that the space density of bright 
red galaxies is nearly constant over 1.5 > z > 0.5, 
while that of bright blue galaxies decreases by a 
factor of ~ 2 over the same redshift interval. 

Recently, lllbert et al.1 (l2010h found a rise in the 
stellar mass density of log M^/Mq > 11 galaxies 
from z~2toz~lbya factor ~ 14 for quies- 
cent galaxies and a factor of ~ 4.3 for red-sequence 
galaxies, similar to our result. This difference be- 
tween quiescent and red-sequence galaxies is likely 
to arise from red star-forming galaxies that con- 
taminate the red sequence more towards z ~ 2. 
For both quiescent and red-sequence galaxies, they 
found little evolution at z < 1, as we do. Also, 
their highly star-forming sample compares well to 
our blue galaxy population. They find a rise in 



mass density by a factor of ~ 1.5 from z ~ 2 to 
z ~ 1 and a decrease by a factor of ~ 2.5 at lower 
redshift again. 

6. SUMMARY AND CONCLUSION 

In this work we investigated the evolution of 
the red sequence in terms of colour, luminosity, 
mass, and number density through cosmic time 
since z = 2. We derived an ii-band catalogue of 
10692 galaxies from 0.2 deg 2 of the A901-field sur- 
veyed by the deep NIR multi-wavelength survey 
COMBO-17+4. While deep multi-wavelength sur- 
veys provide photometric redshifts for large sam- 
ples of galaxies, the measured colours suffer from 
large uncertainties at high redshift that wash out 
the contrast with which the red sequence appears 
on top of the tail extending from the blue cloud. 
Hence, we used colour histograms weighted by the 
colour error of each galaxy to trace even a diluted 
red-sequence signal. We also used the rest-frame 
colour (C/280 — V) for maximum population con- 
trast and minimum sensitivity to redshift errors. 

As a result, we found a red sequence up to 
z ~ 1.65, beyond which the situation is unclear. 
Tracking the colour evolution of the red sequence 
peak, we derived an evolving colour-magnitude re- 
lation up to z — 2 and used it to separate the red 
and blue galaxy populations. Our results show 
that the (l^so — V) colour evolution of the red 
sequence is consistent with pure aging. Further 
results are: 

1 . Both the red and blue galaxy population get 
redder by A(U2so — V) ^0.4 mag since 2 = 2. 

2. The population of massive blue galaxies 
grows from z ~ 2 to z ~ 1. 

3. The massive end of CM*Ds is dominated by 
the red galaxy population at z < 1 and by 
both galaxy populations at z > 1. 

4. Some massive red galaxies with log M* /M & ~ 
11.5 are already in place at z ~ 2. 

We investigated the number density and stellar 
mass density evolution of massive red galaxies and 
found that both increase by a factor ~ 4 between 
2 > z > 1 and shows little evolution since z = 1. 
This suggests that the main formation epoch of 
massive red galaxies is at 2 > z > 1 such that 
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they have assembled most of their mass by z ~ 1 . 
Note, that our masses are unlikely to be biased 
much by dust. 

It is clear that our results are affected by cosmic 
variance due to the small area surveyed. Once the 
data are available from the two other fields (A226 
and Sll) targeted by the COMBO-17+4 survey, 
we expect to create an if -band selected catalogue 
of ~ 50000 galaxies in an area of 0.7 deg 2 , of which 
~ 12000 galaxies will be above z = 1. The more 
than threefold increase in area and the combina- 
tion of unrelated fields on the sky will reduce the 
effect of cosmic variance by a factor three and firm 
up our findings quantitatively. 

We are grateful to all Calar Alto staff as- 
tronomers in particular Jesus Aceituno, Anna 
Guijaro, Felipe Hoyo and Ulli Thiele. Thanks 
to the International Max-Planck Research School 
for fundings. CW was supported by an STFC 
Advanced Fellowship. 
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Fig. 1. — Spectral energy distributions of four red- 
sequence galaxies at different redshift. 
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Fig. 2. — Photometric redshift accuracy vs. H- 
band magnitude in higher-z redshift intervals. 
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Fig. 3. — Photometric redshift distribution of our 
galaxy sample vs. H-band magnitude. 



13 




Fig. 4.— Left: Original COMBO- 17 photo-z dis- 
tribution in four fields, illustrating field-to-field 
variation. Right: New COMBO-17+4 distribution 
of the sample used in this paper (black line), and 
the red sample as defined in Sect. 4.3 (grey line). 
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Fig. 5. — Panel a): Rest-frame colour (7/280 — V) 
over redshift. The A901/2 super- cluster produces 
a prominent feature at z—0. 165. Panel b): Colour 
error over redshift, as determined from error prop- 
agation. At high redshift the colour determination 
is less accurate, especially at z > 1.5. An error 
floor o/0™l is assumed, see Sect. 
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Fig. 6. — Rest-frame (C/28O — ^)m v =-20 galax- 
ies colour distribution in redshift slices. His- 
tograms are sums of Gaussian probability distri- 
butions representing single galaxies with their in- 
dividual colour error. Structure is driven by bright 
objects with accurate colours while faint, large- 
error objects with broad, smooth distributions don't 
add contrast to the peaks. The top left panel con- 
tains the super-cluster A901/2 with its prominent 
red sequence. The colour of the red sequence peak, 
indicated by an arrow, is determined by a Gaus- 
sian fit to the localised excess above a continuum 
of galaxy counts. Thus, a red sequence is observed 
up to z ~ 1.6. Number of galaxies, mean redshift 
and mean lookback time t — z/(l + z) are noted. 
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Fig. 7. — Evolution of the red-sequence colour with 
lookback time t — 15 Gyr ■ z(l + z)^ 1 , see Table 
[H The dashed line is a linear fit to the points, 
and the solid line is an example prediction using 
PEG AS E for a single age stellar population with 
solar metallicity formed 12 Gyr ago (zf — 3.7). 
The error bars account only for the colour mea- 
surement of the peak of the bright end of the red 
sequence. However, the scatter among the points 
illustrates further uncertainties including system- 
atics and cosmic variance. 



15 




0.94<Z<1.07 Total=1077: 
Rcd=266 



1.07<Z<1.19 




; 1.55<Z<1.78 
3 — Volume *2 5 



; 1.78<z<2.0 
3 - Volume x25 




Fig. 8. — Colour-magnitude diagrams in redshift 
slices. Solid lines indicate the bimodality sepa- 
ration (see Eq. [5)j at the mean redshift of each 
bin. The real separation is given by the colour of 
the data points (black=red sequence, green=blue 
cloud). Total number of galaxies and red fraction 
are noted. All redshift bins are nearly equal co- 
moving volumes except for the three highest bins, 
which enclose 2.5 times as much volume, and the 
A901/2 cluster bin with a small volume enclosed 
within 0.150 < z < 0.175. For reference, we show 
a reddening vector for Ay = 1 mag. (A colour 
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Fig. 9. — Colour-mass diagrams in same redshift 
slices and for mat as Fig. O The re ddening vec- 



tor assumes a \Bell & deJonal \200A) law fo r the 
stellar-mass- colour relation and a \Pei I 19921) dust 
law. (A colour version of this figure is available 
online.) 




Redshift 

Fig. 10. — Number density evolution of massive 
red galaxies with log M»/Mq > 11 and (E/28O — 
V) > (U2S0 ~ V)um as a function of redshift. The 
error bars represent the cosmic variance in each 
redshift bi n and were calc u lated with a method de- 
scrib ed bu \Moster et al. I If 20 id) . 



from \Tavlor et all \200d) /< 



In 



grey, 



data 



for comparison. 
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0.5 1 1.5 2 

Redshift 

Fig. 11. — Stellar mass density evolution of galax- 
ies with log M*/M© > 11 as a function of red- 
shift. The entire galaxy population (diamond) is 
divided into the red-sequence (square) and blue- 
cloud galaxy population (triangle). Error bars in- 
clude only cosmic variance estimates and no un- 
certainty from the stellar mass estimation. (A 
colour version of this figure is available online.) 
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Table 1 

COMBO-17+4 Target Coordinates and Integration Times 



Field a 


RA 


DEC 


Integration 


Time per 


Filter 13 








H 


J 2 


Jl 


Y 




(J2000) 


(.T2000) 


(ksec) 


(ksec) 


(ksec) 


(ksec) 


A901 


09h 56m 17s 


-10° 01'25" 


11.6 


11.9 


8.6 


8.4 


A226 


Olh 39m 00s 


-10° n'oo" 


14.6 


17.3 


11.6 


8.8 


Sll 


llh 42m 58s 


-01° 42'50" 


16 


4.5 


8 


10 



a Rcsults in this paper are derived from the observations of the A901-ficld only, 
integration time averaged over different subfields. 



Table 2 

(£^280 - V)m v =-2o Colour of the Red-Sequence Peak 



Mean Rcdshift 


Redshift Interval 


(U28O — V)m v = -20 

(Vega mag) 


<5 (^280 — U)m v = -20 

(Vega mag) 


Lookback Time 
(Gyr) 


0.165 


0.15-0.175 


2.10 


0.04 


2.034 


0.358 


0.3-0.4 


1.704 


0.03 


3.923 


0.525 


0.45-0.6 


1.70 


0.03 


5.217 


0.720 


0.65-0.78 


1.48 


0.03 


6.432 


0.821 


0.78-0.86 


1.27 


0.03 


6.960 


0.923 


0.86-0.97 


1.045 


0.03 


7.437 


1.097 


1.05-1.15 


1.1 


0.03 


8.139 


1.246 


1.2-1.3 


1.05 


0.03 


8.647 


1.397 


1.3-1.5 


0.68 


0.03 


9.091 


1.552 


1.5-1.6 


0.85 


0.03 


9.486 


1.646 


1.6-1.7 


0.59 


0.03 


9.700 



Table 3 

Number and Stellar Mass Density of Galaxies with log M*/M q > 11 



Rcdshift Interval 


4>red 






lo S P'all 




(10~ 4 Mpc~ 3 ) 


(M Q Mpc" 3 ) 


(M Q Mpc~ 3 ) 


(M Mpc~ 3 ) 


0.2-0.6 


2.550 


7.627 


6.825 


7.691 


0.6-0.79 


2.911 


7.671 


6.563 


7.704 


0.79-0.94 


3.169 


7.765 


6.932 


7.824 


0.94-1.07 


1.890 


7.508 


7.117 


7.656 


1.07-1.19 


2.555 


7.648 


7.167 


7.825 


1.19-1.3 


1.984 


7.538 


7.241 


7.716 


1.3-1.55 


1.323 


7.438 


7.025 


7.579 


1.55-1.78 


1.125 


7.337 


7.229 


7.586 


1.78-2.0 


0.643 


7.112 


6.818 


7.346 
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